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Abstract 

We get the leading term of the Gurevich-Pitaevskii special solution of the 
KdV equation in the oscillation zone without using averaging methods. 

Keywords: nondissipative shock wave, cusp catastrophe 
1. Introduction 

The Gurevich-Pitaevskii (GP) special universal solution of the Korteweg- 
de Vries (KdV) equation 

ih + u u x + u xxx = (1.1) 

was introduced in [l[ in connection with the problem of description of colli- 
sionless shock waves (Sagdeev showed in [2J that such waves are of oscillat- 
ing character). The behaviour of the GP special solution for t — > — oo and 
x — > ±oo is determined in the main from the cubic canonical equation of the 
cusp catastrophe 

x-tu + u 3 = 0. (1.2) 

The GP solution to the KdV equation is one of the most interesting 
special functions of the modern nonlinear mathematical physics. 
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In [l| it is shown that in problems of dispersion hydrodynamics (in par- 
ticular, in problems of plasma theory) the GP special solution appears near 
the points of overturning of simple waves. From the results of Q, jij, 
one actually sees that the same universal special function appears near the 
points of overturning of the generic state solutions to diverse dispersion per- 
turbations of the equations of one-dimensional motion of ideal incompressible 
liquid 

Pt + (vp)' x = 0, 
v' t + v v' x + a(p) p' x = 0. 

Here, p is the density of the liquid, v the velocity and a(p) = (c(p)) 2 /p, where 
c (p) = yP'(p) is the speed of sound and p(p) the pressure. In particular, 
this is the case for solutions of the shallow water equations 



K + (hA' x )' x = e 2 (hZA» x )» x + 0(e*) 
\(A' x f + gh = \ 



A' t + -(A' x ) 2 + gh = -e 2 (A:> xt + A> x A:> xx -(A: x r) + 0(e 4 



where h is the free boundary, A the potential of bottom velocity and g 
the acceleration of gravity. The right-hand sides can actually be written 
as complete series in powers of the parameter e by the procedure given for 
instance in [6, Ch. 1, §4] (and not only as the so-called second approximations, 
as stated in 4j). 

In the 1990s there were discovered surprising connections of the GP spe- 
cial solutions with some problems of quantum gravity. In [t| this solution 
was showed to simultaneously satisfy the fourth order ordinary differential 
equation 

5 5 5 

u xxxx + -uu xx + -(u x ) 2 + — (x-tu + u 3 ) = 0, (1.3) 

3 o 18 

which had been studied for t = in js|] and ^ in connection with evaluat- 
ing nonperturbative string effects in two-dimensional quantum gravity (the 
equation (11. 3ft belongs to a class of massive string equations). In [10] the 
solution of 

(W 3 - WW XX - ^{W x f + ^Wxxxx^j + ^T(W 2 - Uv X x) = X 

with asymptotics \fX as X — > ±oo was treated numerically in connection 
with problems of quantum gravity. One can show that this solution U(t, X) 
is also equivalent to the GP special solution of (II. ip for t > (but not for 
t < 0). 
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Dubrovin showed in and 12[ directly by means of the theory of ap- 
proximate symmetries [13J that it is the solution of (II. 3p with asymptotics 
(11.2j) that appears near the points of wave overturning for the very diverse 
singular dispersion perturbations of the equations of one-dimensional hydro- 
dynamics. 

The results of numerical simulations presented in 10] demonstrate rather 
strikingly that the GP special solution of the KdV equation possesses a do- 
main of undamped oscillations for t large enough. The authors of [lOj did not 
conjecture any relation of their paper to the GP special solution and raised 
the problem of describing this domain of oscillations. Meanwhile Gurevich 
and Pitaevskii [14j had used successfully the self-similar solutions of the av- 
eraged Whitham equations 15| to solve the problem. 

The self-similar solutions in question were constructed in explicit form 
by Potemin jig . However, the problem on the leading term of asymptotics 
of the GP special solution in the domain of Whitham oscillations has been 
open up to now. One not simple question still unanswered has been that on 
the phase shift. 

Our purpose is to show how it is possible to construct the leading term of 
the GP special solution in the zone of oscillations without using any averaging 
methods. To this end we derive certain algebraic equations for the slowly 
varying amplitude and the leading term of the phase, which are actually 
equivalent to those of [16]. Moreover, we determine the phase shift of the 
solution in the oscillation zone. 

Our approach may also be of use for the study of undamped oscillations of 
other common solutions to integrable partial and ordinary differential equa- 
tions which are of importance in physics. In particular, it applies to two 



universal solutions of the KdV equation treated in the recent article [17 
Almost one problem in the approach is some awkwardness of analytical cal- 
culations. However, invoking modern programs for symbol calculations (in 
this paper we use Maple) often allows one to get rid of such problems without 
particular difficulties. 



2. Evaluation of phase shift 

Consider the solution of the KdV equation that, for t — > — oo and x — > ±oo, 
is determined in the main from the cubic equation (\1.2\i . It is known that 
for this solution for positive t there is a domain where dissipationless shock 
waves appear. 
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We are aimed at constructing asymptotics of the solution in this domain, 
when t — > oo. Following familiar techniques, we change the variables by 



u = ^/\t\U(t,z), 

X 

Then equations (II. ip and (11. 3p take the form 



tu t + -{u - 3zU z ) + uu z + r 7/2 u z 



0. 



2 5 " (2.1) 

t- 7 U zzzz + ^-t- 7 / 2 (6UU zz + 3(U z ) 2 ) + ^-(z-U + U 3 ) = 0. 

lo lo 

We now look for a solution U of the system in the form of asymptotic 
series 

U = U {cp, z) + t-'^U^, z) + t- 7 / 2 U 2 (p, z) + ..., (2.2) 

where Uq, U\ and U 2 are 2ir -periodic in the fast variable ip. This latter is 
assumed to be of the form 

p = t- 7 /*f(z)+s(z), 

where by s(z) is meant precisely the phase shift. 

For the unknown function Uq we get the nonlinear system 



Q 3 d 3 U Q + QRd v U + QU d v U = 0, 
Q^U + ^Q 2 (2U d 2 v U + (d v Uo) 2 ) + ^(z-U + U 3 ) = 0, 

while the systems for U% 

Q^Ui + Q(R + U )3 V U 1 + Qd^Uo U x = F u 
Q%U X + h - Q\V^%V X + dyUadyUi) + ^ (3U 2 + 3Q 2 dJf/ - l)U x = F 2 , 

and for U 2 



Q 3 d 3 U 2 + Q(R + U )d v U 2 + Qd^Uo U 2 = G u 
^Q 2 (U dlU 2 + d ip U d v U 2 ) + ^ 



Q A d^U 2 + - Q 2 (U Q dlU 2 + d v U d v U 2 ) + — (3U 2 + 3Q 2 d 2 v U - l)U 2 = G 2 



proves to be linear. Here, F\ and F 2 are explicit functions depending on z 
and Uq, and Gi and G 2 are explicit functions depending on z and U , U\, 
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i.e., the right-hand sides are explicit functions depending on z and on the 
preceding corrections. We write 



Q = f, 

7 

4f 2 



R = U- 3 - Z (2 ' 3) 



for short. 

From the compatibility condition of the equations for Uq we obtain a first 
order equation 

Q 2 (d^U ) 2 + + R U 2 + \ (18R 2 -h)U Q + \{lhR- 54i? 3 - 5z) = 0. (2.4) 
o o o 

From the compatibility condition of the equations for U\ we derive a nonlinear 
equation 

d 1 486 R 4 - 171 i? 2 + 9zR + 5 

di ~ 9 (54 J R 3 -9 J R + ^)(2 J R + 3z) ^ ' ' 

for the unknown function R = R(z). (In Section H] we show that this equation 
agrees with results obtained earlier.) When requiring the compatibility of the 
equations for U 2 , we deduce that the function U (ip, z) should satisfy, together 
with (12.41) . a nonlinear ordinary differential equation in the variable z of the 
form 

^•-(w (w+^-w. (s " +ffs ' )H w = °- (2 - 6) 

Here, Ps(U ) and P±(Uq) are polynomials in U of degrees 3 and 4, respec- 
tively, with coefficients depending on z and R. The function H = H(z, R) is 
given by 

h(z m = 1 - N ^ R) 

1 ' ; 3 (54R 3 -9R + z) 2 (2R + 3z) 2 

where 

N(z,R) = 45z 3 + (4860R 3 -582R)z 2 + (131220 J R 6 -43416 J R 4 + 2721i? 2 -35)z 
+ 139968 J R 7 -59616 J R 5 + 6048 J R 3 -120 J R. 

Note that (I2.6P is a Hamiltonian equation with Hamiltonian quadratic 
relative to the impulse, i.e. a(d z Uo) 2 + bd z Uo + c where a, b and c are functions 
of z and Uq. 
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We proceed to study equations fl2.4p - fl2.6p . Equation (12. 4p is autonomous 
in the fast variable ip, hence the arbitrary constant of the general solution is 
contained in the phase shift which we take into account in the variable s(z). 
The general solution of this equation is sought in the form 

U = Adn 2 ^k)+C, 

where dn is the elliptic function of Jacobi and A, B, C and k are to be de- 
fined. On substituting Uq into equation (12 .4p and equating the coefficients of 
different powers of the Jacobi function to zero we get the system of algebraic 
equations 

A 2 {A- 12B 2 ) 
A 2 ((8-2k 2 )B 2 + C + R) 
A ({12k 2 - 12)AB 2 + 18A 2 + 6RC - 5) 
C 3 + 3RC 2 + (18R 2 - 5)C - 54R 3 + IhR - hz 

From the assumption on the 2n -periodicity of Uo it follows that 

| = Z®, (2.8) 

where K{k) is the complete elliptic integral of first kind. 

The system of equations (12. 3ft . f 12 .51) . ( 12. 7ft and ( 12.81) obtained in this 
way is overdetermined. It consists of 6 algebraic equations and 2 differential 
equations for the unknowns /, Q, R, A, B, C and k. Our next concern will 
be to show that one can find all slowly varying unknowns without solving 
the differential equations. 

The unknowns /, Q, R, B and C can be determined immediately from 
this system through A, k and z. More precisely, 

(k 2 -2)(4k A -5k 2 + 5)A 3 -l0(K(k)-2)A-4:5z 
14(fc 4 - k 2 + l)A 2 - 30 ' 

TlB 

K{kj' 



(2.9) 



0. 

0. 
0. 
0. 



(2.7) 



B 
R 

f 




{k 2 -2)A 
3 

Q(AR+6z) 
7 



C 



■C, Q 
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On substituting these expressions into (12.71) we arrive at one algebraic equa- 
tion 



(8A: 12 - 24k 10 +43 A; 8 -46 k 6 - 43 k 4 + 24k 2 - 8) A 6 - 140( k 4 - k 2 + l) 2 A 4 
+ Wz(k 2 -2){2k 2 -l)(k 2 + l) A 3 + 500 (A; 4 - A; 2 + 1) A 2 +3375z 2 = 500 

(2.10) 

for A and k. Differentiating this equality in z and substituting the result- 
ing expression along with (12 .9p into (12. 3p and (12. 5p . we obtain 3 equations 
containing A'(z) and k'(z). On eliminating these derivatives, we get another 
algebraic equation for A and k which contains the quotient of two complete 
elliptic integrals 

q K{k) 

Using (12.101) to eliminate the highest power of z from the latter equation we 
bring it to the form 

21A: 4 (A; 2 -l) 2 v4 3 + 10((2g-l)A; 6 -(3g+l)A; 4 -(3g-4)A; 2 + (2g-2))A 
+ 315,2 ((2g-l)A; 4 -(2g-3)£; 2 + (2g-2)) = 0. 

(2.11) 

When eliminating the variable A from (I2.10p and (12. lip , one finds k as an im- 
plicit function of z. The other functions can be expressed explicitly through 
k. 



This method allows one to get the explicit formulas of [16[ without using 
averaging procedure. 

The domain of z in which oscillations are possible is determined from the 
condition that k G [0, 1]. To the point k = there corresponds leading wave 
front set zi = —a/2 and to the point k = 1 there corresponds trailing wave 
front set z t = VTO/27. 

The second ordinary differential equation of (12.11) enables us also to de- 
termine the phase shift s(z) of the solution. For this purpose we make use 
of equation (12. 6p . Note that Uq is an even function of (p, hence d^Uo is odd 
and d z U , d 2 U are even in <p. Multiplying (12. 6p by (d v U ) 3 , integrating in 
over the whole period 2n and taking into account that the mean value of an 
odd periodic function over the whole period vanishes, we get the equation 

s" + H(z,R)s' = 0. (2.12) 

In order to choose a concrete solution to (I2.12p . one has to use the asymp- 
totics of the solution for z — >■ z\. Such an asymptotics is given in [i[ and it 
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shows in particular that the solution does not contain log(z-zi). In our case 
we get 

. 2 7 /8 ^ , 2 V\ t1 „ 131 



V5 



(^/) 1/4 (l-^(^) 1/2 + ^V2(^)+0((^) 3/2 ))^ 



2—2;; 1600 

From the asymptotics of H it follows that a solution of ( I2.12p contains 
terms of the form log(z + v^2) which are not permitted by jj Hence, this 
solution enters into the linear combination of solutions with coefficient 0, and 
so 

s(z) = s 

is constant. 

The hypothesis on the constancy of the phase shift s(z) has been formu- 
lated in [if. However, in (jjj it was based solely on the asymptotics given 
there. It is clear that a priori one might not exclude the situation where the 
phase shift fails to be constant but tends exponentially fast to a constant as 
z — > Z\. From (I2.12p and the asymptotics of H(z, R) for z — > Zi we see that 
such is not the case, and so s(z) proves to be constant. 

To evaluate the constant sq we invoke a numerical simulation. Namely, 
we compare a numerical solution with the solution constructed by using 
asymptotic formulas. 



3. Numerical simulations 

To this end we have written a special program. The results of numerical 
simulations are presented in Figures [1] - [3j 

In Figures [T]|2] one can observe numerical solutions for function U(t, z) 
for negative and positive value t. It can be shown that function U(t, z) for 
negative value t practically councide with the root of the equation A 3 + A = 
—z. 

In Figure [3] the difference between these solutions is shown. Both figures 
correspond to t — 20. The constant so has proved to be equal to 3.1254 w n. 
The difference between two solutions is a multiple of t -L77 w i~ 7 / 4 which 
agrees with the order of the first correction in formula (12. 2p . 
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Figure 1: The numerical simulation for the function U(t, z) corresponding to t = — 7. 



We believe that the constant Sq just amounts to ir and the small difference 
is caused by a computation error. 

Figure [3] makes it evident that the error increases for z close to — v2 and 
VTO/27. This manifests the nonuniform character of the constructed asymp- 
totic formula in the entire domain of the variable z. In neighbourhoods of the 
leading and trailing wave front sets one should construct other asymptotic 
formulas which can be made consistent with asymptotic expansion (12. 2p . as 



it is described in 18 



4. Reduction of equations to standard form 

To start with, we bring the equations obtained in |l6j to a simpler form. 
These are 

z-Yj + Zj = 

for j = 1, 2, 3, where 

1 2h-h 1 2 (Z 1 -Z 2 )(l-fc 2 ) 1 2l a -l 2 

Yl ~ 3 1 + 3l^"' Y2 ~ 3 1_ 3 ' Yi ~ 3 1 + 3^~' 

S\ = I1+I2 + h, 5*2 = M2 + hh + W2, S3 = hhh 



and 



with V = 5Sf - 12SA + 8S 3 . 
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Figure 2: The numerical simulation for the function U(t,z) corresponding to t = 20 and 
root of equation z — A + A 3 = 0. 



Obviously, one can eliminate two variables from three equations. On 
eliminating z and q we get an equation for /i, l 2 and l 3 , which does not 
contain z and q. Namely, 

3 (Ij + l 2 2 + l 2 3 )+2 (hl 2 + l 2 l 3 + l 3 h) -5 = 0. (4.1) 

We now eliminate either of the variables z and q in the equations and use 
(14. II) to eliminate all powers of l\ greater than the first one. Then we get two 
more equations 



— (/i(8/2 + Al 2 l 3 + 8% - 15) - (Z 2 + h){m\ - 8/ 2 / 3 + 2Al\ - 25)) . 



2 

45 

1 (l 2 -l 3 )(M 2 l 3 + 3hh + 9ll-5) 

q ~ 2 h(2l 2 2 + l 2 l 3 + 211 - 5) - (k + h)(Ql 2 2 ~ 2hh + 6Z| - 5) ' 

(4.2) 

The system of equations (I2.3p . (12. 7\i and (12. 8p is equivalent to the sys- 
tem of averaged equations (14. ip and (14. 2 p which was obtained in [l6| by 
Whitham's method. To prove this, we pass in the system of the equation for 
C in (J2HD and equations f l210|) . fl2TTT|) to Whitham's variables k(z), l 2 (z) 
and l 3 (z) by 

A = 2(l 3 -h), C = h + k-h, k 2 = ^— l l 

h — n 
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Figure 3: The difference between the numerical solutions and asymptotic for U (t, z) for 
t = 20. 



On eliminating z and g from the obtained equations we arrive precisely at 
(I4.ll) . On eliminating either of z and q and all powers of li greater than the 
first one, we get 04. 2 1) , as desired. 



5. Conclusion 

It should be noted that the value of the phase shift sq derived from the nu- 
merical simulation contradicts 0, 0] • From the results of the paper it follows 
that the function s(z) tends to 7t/2, as z — > Zj. But there is an arithmetical 
error m formula (17) (jj and in formula (21) (J). The right calculation with 
using monodromic date from [19J showed that s(z) tends to ir, and so 

s(z) = it. 
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